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Abstract 

The problem of factorization of a nonextensive probability distribution is discussed. 
It is shown that the correlation energy between the correlated subsystems in the 
canonical composite system can not be neglected even in the thermodynamic limit. 
In consequence, the factorization approximation should be employed carefully ac- 
cording to different systems. It is also shown that the zeroth law of thermodynamics 
can be established within the framework of the Incomplete Statistical Mechanics 
(ISM). 
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1 Introduction 



The nonextensive probability distribution 

„ _ [i-{i-q)m-c)\^ 

plays a decisive role for the success of the generalized statistical mechanics [1,2] 
because it is capable of reproducing unusual distributions of non gaussian type 
which are met frequently in nature (see reference [1] and the references there- 
in). Very recently, distribution functions of type Eq. (1) were recognized by 
the so-called eigencoordinates method with high level of authenticity as good 
description of the amplitude distribution of earthquake noises which are proved 
to be fractal and strongly correlated [3]. In Eq. (1), /3 is the generalized inverse 
temperature, Ei is the energy of the system in the state i, C a constant to 
assure the invariance of the distribution through uniform translation of energy 
spectrum E^ and Z q is given by 

Z q = J2ll-(l-q)m-U q )]^ (2) 



or by 

^ = £[l-(l-s)/«-^)]^,[4] (3) 

i 

in Tsallis' multi-fractal inspired scenario [1] with U q the internal energy of the 
system and (3 = and by 

i 

9 

(4) 



£[1 - (1 - <z)/H^ 



with (3 = ^jt- in Wang's incomplete statistics scenario [2] devoted to describe 
inexact or incomplete probability distribution [5] due to neglected interac- 
tions in the system hamiltonian. It is noteworthy that Eq. (1) is a canonical 
distribution function for isolated system in terms of its total energy Ei. 

In this paper, we discuss the factorization problem of Eq. (1) and some of 
its consequences when the canonical system is composed of correlated sub- 
systems or particles. At the same time, we also comment on some interesting 
applications of Eq. (1) found in the literature. 
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2 Factorization approximation 



In Boltzmann-Gibbs-Shannon statistics (BGS), thanks to the easy factoriza- 
tion of the exponential distribution functions, i.e., 

Ee-^^nEe-**. (5) 

i j i 

the system distribution function can yield one particle one in terms of the 
single body energy in the case of independent particles or of mean-field 
method with E^ = J2f=i e iji where j is the index and N the total number of 
the distinguishable particles. 

But this approach is impossible in the case of the nonextensive distribution, 

because, 

Et 1 " (1 " ff)0Eey]^ ^ I1E[1 " (1 " Q)fa]^- (6) 

i j j i 

This inequality makes it difficult to apply the nonextensive distribution even 
to systems of independent particles, because the total partition function can 
not be factorized into single particle one. So in the definition of the total sys- 
tem entropy S q = f[p(Ei), q], p(Ei) must be the probability of a microstate of 
the system and can not be replaced by one-body energy e^. However, in 
the literature, we find applications of Tsallis' distribution in which Ei is sys- 
tematically replaced by ejj without explanation [6-11]. The first examples[6,8] 
are related to the polytropic model of galaxies and the authors have taken 
the one-body energy of stars and of solar neutrinos (e = ^ + y ) as Ei. Other 
examples are the peculiar velocity of galaxy clusters (e ~ ^ 2 )[H] and the elec- 
tron plasma turbulence where the electron single site density n(r) was taken 
as system (electron plasma) distribution function and the total energy was 
calculated with the one-electron potential 4>(r) [7,9]. The last case is the ap- 
plication of nonextensive blackbody distribution to laser systems where the 
atomic energy levels were taken as laser system energy in equation (8) of ref- 
erence [10]. As a matter of fact, in above examples, the calculated entropies 
and distributions are one-particle ones, as we always do in BGS framework. 
Consequently, considering Eq. (6), they are only approximate applications of 
the exact Tsallis' distribution. 

The legitimacy of the above mentioned applications depends on the approx- 
imation with which we write Eq. (6) as an equality. One of the solutions is 
the limit of q — > 1 where Eq. (6) tends to an equality. But this solution does 
not hold for the cases of q value very different from unity, as in the above 
mentioned examples of applications. Another way out is the factorization ap- 
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proximation proposed in order to obtain the nonextensive Fermi-Dirac and 
Bose-Einstein distributions [12] which read 
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(7) 



[l + ( g -l)/3(e-/x)]^±l 



where g is the degeneracy of the level with energy e. As in the case of the 
corrected Boltzmann distribution (i.e. -<1), Eq. (7) can be reduced to 



where p(e) is one-particle probability and e one-particle energy. So we can say 
that only in this approximation the above mentioned applications are justified. 

It is worth mentioning that, although approximate, the above successful appli- 
cations were the first proofs of the existence of Tsallis type distributions (Eq. 
(1)) in nature. In addition, this is a parametrized distribution. So in many 
cases, what we neglect in the approximation may be compensated (at least 
partially) by a different value of q fixed empirically. The interest of the applica- 
tions mentioned above is to show that this kind of nonextensive (nonadditive) 
probability can really describe some non-gaussian type peculiar distributions 
we observe. That is what is important in practice. 

Nevertheless, an approximation has sometimes in itself theoretical importance 
when it concerns the basic foundation of a theory. That is the case of the 
factorization approximation. 



3 Factorization of the joint probability 

The factorization approximation is a forced marriage between the right-hand 
side and the left-hand side of Eq. (6), that is we write just like that : 



where we keep only the index j of the particles, as we do from now on in this 
section. What is neglected in this approximation is the difference A between 
the right-hand side and the left-hand side of Eq. (9) which has been investi- 
gated in reference [13] with a two-level system for the simple case where q > 1 
(or q < 1), ej > (or Bj < 0) and \i = 0. Under these harsh conditions, A 
turns out to be very small at normal temperatures for mesoscopic or macro- 
scopic systems (with important particle number N). But it is not the case for 





[i - (i - ^ IB 1 - a - ?)W^5. 



(9) 



3 3 
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q < 1 and ej > with in addition /i ^ 1. So in general, we can not write Eq. 
(9). 



Eqs. (6) and Eq. (9) can be discussed in another way as follows. If we replace 
J2j ?j by the total energy E in Eq. (9), we get : 

[1 - (1 - q)f3E}^ = W ~ (1 - ?)/fe,-]^5. (10) 

j 

for N subsystems (or particles with energy ej where j = 1, 2, ...TV) of a compos- 
ite system with total energy E at a given state. This is just the factorization 
of the joint probability p(E) as a product of all p(ej) : 

N 

p(E) = Hp(e j ). (11) 

With Eq. (10) or (11), strictly speaking, we can not write 

E = J2e r (12) 

j 

But in Tsallis' scenario, Eq. (12) is necessary for establishing the zeroth law 
of thermodynamics [14]. Abe studied this problem with ideal gas model and 
concluded that Eq. (12) can hold for e,- > 0, < q < 1 and N — > oo. That is 
the correlation energy E c = E — J2j e j between the maybe strongly correlated 
subsystems or particles can be neglected. But this is of course not a general 
conclusion for any q value or any system. In what follows, we will try to give 
the general expression of the correlation energy in ISM because this relation 
is implicit in Tsallis' scenario [15]. The following discussion is for < q < oo, 
the permitted interval of q value in ISM [2] . 

If N — 1, from Eq. (11), we naturally obtain E = e x so E c = 0. 

If iV = 2, we obtain : 

aE = {l + ae 1 ){l + ae 2 ) -I (13) 

2 

= at\ + ae 2 + a 2 eie 2 = a e« + a 2 eie 2 

i=i 

where a = (q — l)/3. So E c = aeie 2 . 
If iV = 3, we get : 
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aE = (l + ae 1 )(l + ae 2 )(l + ae 3 )-l (14) 
= a(ei + e 2 + e 3 ) + a 2 (eie 2 + eie 3 + e 2 e 3 ) + a 3 (eie 2 e 3 ) 

3 3(^\jterms) 3 

= a^e; + a 2 ^ e^e^ + a 3 J| e*. 

and 

£" c = a(eie 2 + eie 3 + e 2 e 3 ) + a 2 (e!e 2 e 3 ). (15) 
When N — > oo, this is a infinite product problem. In general, we obtain : 



TV 

a£=n( 1 + ae - 1 (16) 

i=i 

AT Af ( 2!(JV-2)! termS ) 

i=l «i<«2 

(JV-3)! termS ) W ( iV!(JV-]V)! termS ) AT 

+ a 3 ^e i2 e i3 + ... + a w ^ JJ e< 

U<«2<«3 ii<i'2<...<iN i=l 

=z« fc e ir 

fc=l ii<i2<...<jfc J=l 

= aE'o + ai? c 

where £q is the system energy given by Eq. (12) for independent subsystems 
like the particles of perfect gas. E c is given by : 



N 

E c = H(l + a ei )-E (17) 
i=i 

AT Ar (feT(# L fc)T terms ) fc 



=e«*- 1 e n 

fc=2 ii<i2<...<ifc i=l 



The value of E c is in general difficult to estimate. We can discuss it in the 
following way. 

For q > 1 (or < q < 1) and ej > (or e« < 0), £7 C may be very big when 
N — > oo because X^( ae i) ma Y diverge. So it is in general not negligible. 

For g > 1 (or < q < 1) and < (or > 0), 
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N N ( k\(N-k)> terms ) k 

^=E(-i)^ 1 [(i-g)/9] fc - 1 E Ilk-I (is) 

k=2 ii<i'2<---<ik J=l 

The sign of each term varies alternatively, which may turn out to cancel E c 
when N — > oo. But this is only a possibility. In general, we can not say that 
E c is negligible. 

Considering Eq. (11), the meanvalue of E c is given by following integral in 
phase-space fi(r) : 

U c = j drp(E)E c 

N 

= E(-i)^[(i 

k=2 
N 

= E(- 1 ) fe ~ 1 [(i 

k=2 

The conclusion of this section is that U c is in general not negligible. So one must 
be very careful in using nonextensive distribution for correlated subsystems 
within factorization approximation, especially in the case of the nonextensive 
quantum statistics derived on the basis of this approximation [12]. 



4 Re-establishment of the zeroth law of thermodynamics 

In this section, we will discuss an application of the so called Incomplete Sta- 
tistical Mechanics (ISM) [2], a new version of the nonextensive statistical 
mechanics (NSM) based on the normalization condition J2iPi — 1 where q 
is a positive parameter, q ^ I corresponds to the fact that the probability 
distributions {p^ with incomplete random variables [5] do not sum to one. If 
q — 1, the random variables of the problem become complete and we recover 
the conventional normalization (CN) J2iPi — 1- ISM was proposed as a con- 
sequence of the study of the fundamental theoretical problems [1,4] of Tsallis 
version of NSM with escort probability [1]. The fundamental philosophy of 
ISM is that CN is difficult to be applied to the systems with important com- 
plicated interactions or correlations because, in this case, our "incomplete" 
knowledge about the states of the systems does not permit us to sum all the 
(exact) probabilities. With simpler systems, this human ignorance could be 
neglected and CN holds. 

In our opinion, ISM is not just another form of Tsallis nonextensive statistics, 



(19) 



ki(N-kyt erms ) k 

E II I / 

ii<i 2 <...<ifc j=l 
N ( kHN-k)< terms ) k 

m k ~ l e iwj 

h<i2<---<ik j=l 
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contrary to what some scientists thought. Some of the significant differences is 
discussed in Ref. [4,15]. In what follows, we will re-establish, in a precise way, 
the zeroth law which was established within Tsallis theory only in the case of 
factorization approximation in neglecting the correlation energy. 

In incomplete statistical mechanics, we have following equations for a com- 
posite system containing two correlated subsystems A and B : 

S q (A + B) = S q (A) + S q (B) + ^Z±S q (A)S q (B), (20) 

Eij {A + B) = Ei(A) + Ej(B) + (q - l)/3Ei(A)Ej(B), (21) 
U q (A + B) = U q (A) + U q {B) + (q- l)(3U q {A)U q {B) (22) 

S q = k — J^ 1 + kf3Zf 1 U q , (23) 



and 



Z l-g 

p= -tr- (24) 



In Eq. (20), there is a plus sign "+" before the correlated term. But it was 
minus sign "-" in Tsallis formalism. Eqs. (21) and (22) show the same in- 
dependence of the correlation terms as Eq. (20). They do not exist in Tsallis 
formalism with escort probability [15]. The (/-dependence of the above equa- 
tions is of crucial importance for the establishment of the zeroth law. 

From Eq. (20), a small variation of the total entropy can be written as : 



5S q (A + B) = [1 + ^-±S q (B)]5S q (A) + [1 + l^S q (A)]6S q (B) (25) 



And from Eq. (22), the variation of the total internal energy is given by : 



8U q (A + B) = [1 + ^U q (B)]SU q (A) + [1 + S_lu q (A)]8U q (B) . (26) 

It is supposed that the total system (A + B) is completely isolated. So SU q (A+ 
B) = which leads to : 



8 



SU q {B) 



SU q (A) 



l + ^U q (A)~ l + *±U q (B) 



(27) 



When the composite system (A + B) is in equilibrium, SS g (A + B) — 0. In 
this case, Eqs. (25) and (27) lead us to : 



1 + *±U q (A) dS q (A) _ 1 + *^U q (B) dS q (B) 
1 + q -^S q (A) dU q (A) 1 + *=±S q {B) dU q (B) 



(28) 



With the help of Eqs. (23) and (24), it is straightforward to show that, in 
general : 



1 + azlrj 

= K" 1 (29) 



which recasts Eq. (28) as follows : 



7 i-g ( A \ dS i( A ) 



- HB) 9S q {B) 
1 } dU q {B) 



(30) 



or 



(3(A) = [3(B) (31) 

where (3 is the generalized inverse temperature defined in Eq. (24). Eq. (30) 
or (31) is the generalized zeroth law of thermodynamics which describes the 
thermodynamic relations between different nonextensive systems in thermal 
equilibrium. 



one of the important meanings of Eq. (31) is that the thermal equilibrium of a 
system is now caracterized by (3 but not T. This allows us to have an explicit 
distribution with Eq. (1), which becomes implicit in Tsallis formalism with 
Eq. (2) and (3). 



5 Conclusion 



Our conclusion is that all applications of Tsallis' distribution relative to cor- 
related subsystems or particles should be considered to be valid only under 
the factorization approximation and the condition y < 1. This factoriza- 
tion approximation must be employed carefully because the correlation en- 
ergy between the correlated subsystems is in general not negligible even in the 
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thermodynamic limit (N — > oo). The zeroth law of thermodynamics can be 
established in the framework of Incomplete Statistical Mechanics without any 
approximation. 
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